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Abstract 
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I. INTRODUCTION 


The decays B —* V. v and B —> D^*H v ( t = e,/i,r) play a prominent role in testing 
the Standard Model (SM) and looking for hints of New Physics (NP) in charged-current 
interactions. In the SM scenario a measurement of these decays provides a direct route to 
determine values of the B meson decay constant fs and the semileptonic form factors. They 
also help to determine the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements \V u b\ and 
\V c b\ to a better precision. A puzzling feature of these decays is that there have been some 
recent hints that lepton universality is broken in the tauonic modes of these decays. 

The leptonic and semileptonic modes are difficult to measure experimentally due to the 
presence of a neutrino in the final state. Ideal in this regard are B-factories where a B meson 
pair is generated from the process e + e _ —>■ T(4S') —> BB. One of the B mesons ( B ta , g ) is then 
reconstructed in hadronic or semileptonic modes, while signal decays of the other B meson 
(B s ig) are identified. A new player has entered the game in that the LHCb collaboration has 
been able to identify the semileptonic decays B° —> D*r~u T and B° —> D */UzU in hadronic 


Q. 


0 , 


collisions 

Since the first evidence reported by Belle collaboration in 2006 [2|J, many measurements 
of the branching fraction B(B~ —> t~u t ) have been reported by both Belle and BABAR 
collaborations. There had been a consistent excess compared to the SM prediction until 
Belle published their result of B(B~ —> t~u t ) = [(7.2^|g(stat) ± l.l(syst))] x 10~ 5 with 
a significance of 3.0cr jsj. This result reduced the tension between theory and experiment 
and decreased the world average of the measured branching fraction to the recent value of 
B(B~ —>• r~u T ) = (11.4 ± 2.2) x 10 -5 ( 4 ], which is slightly larger than the SM expectation 
(8.1 ± 0.7) x 10 -5 obtained from a global fit to CKM matrix elements jj]. Note that the 
most recent result of B{B~ —>• t~u t ) = [12.5 ± 2.8(stat) ± 2.7(syst)] x 10 -5 B reported by 
Belle in September 2014 is in good agreement with its previous result. 

The SM calculation of the leptonic decays suffers from uncertainties in the input values 
of /b and V u b. One can eliminate the V u b dependence by calculating the ratio of branching 
fractions 

_ Tb q B(B~ -» t~u t ) 

w t b - B{B° Tr+£-v e ) ’ 1 J 

where i = n, e. The ratio is measured to be (0.73 ± 0.15) [6|, which exceeds the SM 
prediction of R7 = 0.31 ± 0.06 [6] by more than a factor of 2, while the measured value of 
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B(B° —> 7 = (14.6 ± 0.7) x 10 -5 [7 9 ] is consistent with the SM expectation. 

The semilcptonic decays B —> D'*>lv have a much richer structure than the leptonic 
decays. There is a large number of observables in these decays, e.g., the forward-backward 
asymmetry of the charged lepton. Recently there has been much interest in the ratios of 
branching fractions 

RfRO nWr-Dh 

( 2 ) 




B(B° ->■ D^e-9i) ' 

In taking these ratios some of the uncertainties in the form factors are reduced. Furthermore, 
the dependence on the poorly known CKM matrix element \V c i\ drops out in the ratio. 
Recently, three groups have reported measurements of these ratios 


R(D)\ babar = 0.440 ±0.072 
R(£>)|belle = 0.375 ±0.069 


R(D* 

R(D* 

R(D* 


Ibabar = 0.332 ± 0.030 [10] 

(belle = 0.293 ±0.041 [11] 

|LHCb 


= 0.336 ±0.040 [1] 


where the statistical and systematic uncertainties have been combined in quadrature. These 


measurements were combined in 


12 ] 


R(D) | expt = 0.388 ± 0.047, R(D*) | expt = 0.321 ± 0.021, 

and compared with the SM expectations given in 


(3) 


R(D) Ism = 0.297 ± 0.017, R(D*) | SM = 0.252 ± 0.003. 


(4) 


It is seen that there is a discrepancy of 1.8 a for R[D) and 3.3 a for R(D*). 

The deviation of leptonic and semileptonic tauonic B meson decays from SM expectations 
has been the motivation of many theoretical studies in search for NP effects, including the 
two-Higgs-doublet models (2HDMs) 16-19], the minimal supersymmetric standard model 


(MSSM) 20], and lcptoquark models 


21 


221 ]. In many studies, a general effective Lagrangian 


for the b —>- ubv and the b —> c£v transitions in the presence of NP is imposed to investigate 


li|23|-l25] 


various NP operators and their coupling, together with their correlations 

In this paper we focus on these decays within the SM framework using results from 
our covariant constituent quark model for the dynamics of the transitions. Most of the 
theoretical studies on the semileptonic decays have been relying on elements of the Heavy 
Quark Effective Theory (HQET) 26, ]27|, based on a systematic 1/mQ-expansion of the 
QCD Lagrangian. The leading order of the HQET-expansion corresponds to the Heavy 
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Quark Symmetry when the heavy quark mass tends to infinity, simplifying the structure 
of the weak current transitions. The form factors of these transitions are then expressed 
through only a few universal functions. Unfortunately, HQET can give predictions only 
for the normalization of the form factors at zero recoil. As one moves away from the zero- 
recoil point one has to take recourse to full nonperturbative calculations. In this paper, we 
present a description of these decays that does not rely on HQET. We employ the covariant 
constituent quark model (CQM) with built-in infrared confinement which has been developed 
in several previous papers by our group (see 28j, [29] and references therein). In the CQM 
approach, the entire physical range of momentum transfer is accessible. This is one of those 
features that make the CQM different from other model approaches for the calculation of 


hadronic quantities. We mention that a similar study was done by authors of 


30 


321 


m 


the framework of a relativistic quark model based on the quasipotential approach, in which 
the full range of momentum transfer is also achievable. Our aim is to give an independent 
calculation of these decays including the q 2 behavior of the transition form factors, the 
leptonic decay constants of the B and D mesons, the forward-backward asymmetry of the 
lepton and other polarization observables as well as ratios of branching fractions. 


II. MODEL 


The CQM is based on an effective Lagrangian describing the coupling of a hadron H to 
its constituent quarks, the coupling strength of which is determined by the compositeness 
condition Zh = 0 33j, 34], where is the wave function renormalization constant of 
the hadron H. Here is the matrix element between a physical particle state and the 
corresponding bare state. For ^ = 0 it then follows that the physical state does not 
contain the bare one and is therefore described as a bound state. This does not mean that 
we can solve the QCD bound state equations but we are able to show that the compositeness 
condition provides an effective and self-consistent way to describe the coupling of a particle 
to its constituents. 

One starts with an effective Lagrangian written down in terms of quark and hadron 
variables 3a, [361]. Then, by using Feynman rules, the S'-matrix elements describing hadronic 
interactions are derived from a set of quark diagrams. In particular, the compositeness 
condition enables one to avoid a double counting of hadronic degrees of freedom. This 
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approach is self-consistent and all calculations of physical observables are straightforward. 
There is a small set of model parameters: the constituent quark masses, the scale parameters 
that define the size of the constituent quarks distribution inside a given hadron, and the 
infrared cutoff parameter A. 

The coupling of a meson M to its constituent quarks q\ and q 2 is given by the Lagrangian 

Ant (t) = 9mM(x ) • J M (x) + H.c., (5) 

where cjm denotes the coupling strength of the meson with its constituent quarks. The 
interpolating quark current in ([5]) is taken to be 


Jm{x) 



dx 2 F M (x;x 1 ,x 2 ) q2(x 2 )T M qi(xi), 


( 6 ) 


where the Dirac matrix Tm projects onto the relevant meson state, i.e., Ym — I for a scalar 
meson, Ym = 7 5 for a pseudo-scalar meson, and T m = for a vector meson. The vertex 
function F M is related to the scalar part of the Bethe-Salpeter amplitude and characterizes 
the finite size of the meson. We adopt the following form for the vertex function 


F m (x] xi,x 2 ) = 6 (x- W1X1 - w 2 x 2 )® M ((xi - x 2 ) 2 ), 


(7) 


where wy = rn q J[m qi + m q2 ) so that w\ + w 2 = 1. This form of Fm is invariant under 
the translation Fm(x + a; x± + a, x 2 + a) = Fm(x ; x±, x 2 ), which is a necessary condition to 
provide the Lorentz invariance of the Lagrangian (J5j). 

In order to simplify the calculations we adopt a Gaussian form for the vertex function as 
follows: 

$ M (-p 2 ) = I dxe ipx $ M (x 2 ) = e p2/A M, (8) 


where the parameter A m characterizes the meson size. Calculations of Feynman diagrams 
proceed in the Euclidean region where p 2 = —p%, in which the vertex function has the 
appropriate falloff behavior to provide for the ultraviolet convergence of the loop integral. 

In the evaluation of the quark-loop diagrams we use the free local fermion propagator of 
the constituent quark 


S q {k) = 


m q + ft 


m n 


xe 


mt 


k 2 


te 


with an effective constituent quark mass m q . 


(9) 
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q2 



FIG. 1: One-loop self-energy diagram for a meson. 


For the evaluation of the compositeness condition, we consider the meson mass function 
defined by the diagram in Fig. [0 One has 

n p(p 2 ) = N c g 2 P J + w l p)^S 2 (k - w 2 p)), (10) 

5 v(p 2 ) = N c9v J + w 1 p)j I/ S 2 (k - w 2 p )), ( 11 ) 

where N c — 3 is the number of colors. Since the vector meson is on its mass-shell one has 
ey ■ p = 0 and one needs only the part of the vector meson function proportional to g^ u . It 
is given by 

fivtp 2 ) = - ^)nrw- (i2) 

o p 

The coupling constant gu in Eq. (J5]) is determined by the compositeness condition which is 
written in the form 

Z M = 1 - n' v/ (m 2 v/ ) = 0, (13) 

where II' M (p 2 ) is the derivative of the mass operator taken on the mass-shell p 2 = m 2 M . It is 
convenient to calculate the derivatives of the meson mass functions by using the following 
identities 


o^ n " (p2) - 


W pf W‘ iiMip2) ' 


d 


p M - S(k + wp) = w S(k + wp) $S(k + wp). 

dpv 


(14) 
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Accordingly the derivatives of the meson mass functions can be written as 


x jiui tr [Si(k + wip) i> S 1 (k + wpp)rfS 2 (k - w 2 p)'f 5 ] 
—w 2 tr [Si(k + wip)'y 5 S 2 (k - w 2 p ) piS 2 (k - w 2 p) y 5 ] [, 


( 15 ) 


n W) = 


1 1 


jlV 


p l ‘p v \ 3 g\ f dk ~ 2 


2 p 2 3 y p 2 J 47t 2 ./ 47r 2 i 


(~k 2 ) 


x jwqtr [S^k + Wip) piSi(k + wip)^^S 2 {k - w 2 p)^ u \ 
w 2 tr [S\(k + Wip)^^(k - w 2 p) piS 2 {k - w 2 p) 7„] 


(16) 


The loop integrations in Eqs. (1I5|) and (TT61) are done with the help of the Fock-Schwinger 
representation of quark propagators 

1 rn q + + w pi 


S q (k + wp) = 


rrin 


w pi m 2 — (k + w P y 


= (m q + ft + w pi) Jdae- a[m «-( k+wp)2] . 
o 


(17) 


As will be described later, the use of the Fock-Schwinger representation allows one to do ten¬ 
sor loop integrals in a very efficient way since one can convert loop momenta into derivatives 


of the exponent function (see, e.g., [37H39]). 


As mentioned above, all loop integrations are carried out in Euclidean space. The tran¬ 
sition from Minkowski space to Euclidean space is performed by using the Wick rotation 


k 0 = e*2 ki = iki (18) 

so that k 2 = k$ — k 2 = — k 2 — k 2 = —k E < 0. Simultaneously one has to rotate all external 
momenta, i.e. po — > ipi so that p 2 = —p 2 E < 0. Then the quadratic form in Eq. (fTTli becomes 
positive-definite, 

m 2 ~(k + wp) 2 = m 2 + (k E + wp E ) 2 > 0, 

and the integral over a is absolutely convergent. We will keep the Minkowski notation to 
avoid excessive relabeling. We simply imply that k 2 < 0 and p 2 < 0. 

Collecting the representations of the vertex functions and quark propagators given by 
Eqs. (IS]) and ffTTlh respectively, one can perform the Gaussian integration in the derivatives 
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of the mass functions in Eqs. CB and (fTTTh The exponent has the form ak 2 + 2 kr + zq, 
where r = bp. Using the following properties (k is the loop momentum) 


k M exp(a/r + 2 kr + zq) 
k^k v exp(afc 2 + 2 kr + zo) 

etc. 


1 d 


2 dr. 


exp (ak 2 + 2 kr + zq) 


1 d 1 d 


2 dr M 2 dr u 


exp(aA; 2 + 2 kr + £ 0 ) 


, 


(19) 


one 


can replace ft by $ r — which allows one to exchange the tensor integrations for a 


differentiation of the Gaussian exponent. For example, Eq. (TTOl) now has the form 


n P (p 2 ) = 


3 9 


167r 2 


do,\d(y.2 


tr [ 7 5 (mi + <ft r + wi + <jb r - w 2 ft)\ e ^ +z °. (20) 


o o 


2 / 

The r-dependent Gaussian exponent e~ r ,a can be moved to the left through the differential 
operator <fi r by using the following properties 


d 


dr. 


—r 2 /a g— r 2 /a 


2 d 
+ 


d d , 


dr fi dr v 


V 2 / a = e -r 2 /a 


a drp 

2r fl d 

+ 


a dr.. 


2 r u d 
a dr,. 


etc. 


( 21 ) 


Finally, one has to move the derivatives to the right by using the commutation relation 

d 


dr/ 


= gT- 


( 22 ) 


The last step has been done by using a FORM code which works for any numbers of loops and 
propagators. In the remaining integrals over the Fock-Schwinger parameters 0 < a* < oo 
we introduce an additional integration which converts the set of Fock-Schwinger parameters 
into a simplex. Using the transformation 


OO OO / 

n „ « n « / n 

JJ Idonf(a i,.. ■, a n ) = Idtt n ~ l JJ Idarf M - ^ 
1=1 n n 1=1 \ 1=1 


on ] f(ta 1 ,...,ta n ) (23) 
























one finds 


IW) = 


z 0 

Zm 


OO 1 

0 0 

4 

2s M t 


(24) 


am^ + (1 — a)m q2 — a(l — a)p 2 , 


2 
<?2 

2 2 


(a -w 2 ) P , 


2sm T t 

a M = 2s M +t, b = (a - w 2 )t. 

The function ) arises from the trace evaluation. Further, we have introduced the 

parameter sm = 1/A 2 M . 

ft is readily seen that the integral over t in Eq. (124|) is well defined and convergent if 
Zq > 0, i.e. below the threshold p 2 < (m qi + m q2 ) 2 . The convergence of the integral in the 
case of negative values of zo < 0, i.e. above threshold p 2 > (rn qi + m q2 ) 2 , is guaranteed by 
the addition of a small imaginary to the quark mass, i.e. rn q —> m q — ie , e > 0 in the quark 
propagator Eq. ((2]). It allows one to rotate the integration variable t to the imaginary axis 
t —> it. As a result the integral Eq. (1241) becomes convergent but obtains an imaginary part 
corresponding to quark pair production. 

However, by cutting the scale integration at the upper limit corresponding to the intro¬ 
duction of an infrared cutoff 


oo 1/A 2 

J dt{. ..)—>■ J dt {...), 

o o 

one can remove all possible thresholds present in the initial quark diagram 


(25) 


28]. Thus the 


infrared cutoff parameter A effectively guarantees the confinement of quarks within hadrons. 
This method is quite general and can be used for diagrams with an arbitrary number of 
loops and propagators. In the CQM the infrared cutoff parameter A is taken to be universal 
for all physical processes. 

III. LEPTONIC B-MESON DECAYS 


The model parameters are determined by fitting calculated quantities of basic processes 
to available experimental data or lattice simulations (for details, see Ref. 28], where a 
different set of weak and electromagnetic decays has been used). In this paper we will use 
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the updated least-squares fit performed in Refs. 40|-42]. In this fit we have also updated 
some of the theoretical/experimental input values. The infrared cutoff parameter A of the 
model has been kept fixed. The numerical values of the constituent quark masses and the 
parameter A are given by (all in GeV) 


m n 


m Q 


m r 


m b 


A 


0.241 0.428 1.67 5.04 0.181 


(26) 


Our prime goal is to study the pure leptonic B meson decays as well as the semileptonic 
B —> decays. The most recent results of the fit for those parameters involved in this 


paper are taken from our papers [40 

A d* A d* A d A d 


42|] (all in GeV): 

A Bt A B * 


A b A 


Bs 


A 


B c 


1.53 1.56 1.60 1.75 1.79 1.81 1.96 2.05 2.73 


(27) 


The matrix elements of the leptonic decays are described by the Feynman diagram shown 
in Fig. [2] The leptonic decay constants of the pseudoscalar and vector mesons are defined by 

k + p 



FIG. 2: Quark model diagram for the B meson leptonic decay. 


N c gp 
A J c gv I 


d 4 k 


4 >, 


k 2 ) tr 


(2n)H 

d 4 k ~ . , 9 . 

: $ v (-k 2 ) tr 


(27r) 4 i 


O^S l {k + w l p) 1 b S 2 [k-w 2 p) 
O^Si(k + wip) </.vS 2 (k - w 2 p ) 


fpp 

myfvevi ( 28 ) 


where N c = 3 is the number of colors, and = 7 M (1 — 75 ) is the weak Dirac matrix with 
left chirality. The mesons are taken on their mass-shells. The calculation of the matrix 
elements (I25j) proceeds in a way similar to the case of the mass functions. 

Our results for the leptonic decay constants of Bj*) and D{*J mesons are given in Table [0 
For comparison, we also list the values of these constants obtained from experiments, Lattice 
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and QCD sum rules. Our results show good agreement (within 10%) with results of the 


other studies. We mention that early attempts to account for flavor symmetry breaking in 
pseudoscalar meson decay constants were done in 


43 


44], 



This work 

Other 

Ref. 

I'b 

193.1 

190.6T4.7 

PDG [45] 

I'b s 

238.7 

242.0(9.5) 

LAT [46] 



259(32) 

HPQCD LAT [47] 



193(7) 

LAT [48] 

Ib c 

489.0 

489 ± 4 ± 3 

LAT [49] 

Ib* 

196.0 

196(24)jl| 9 

LAT [50] 



186.4 ±3.2 

QCDSR [51] 

f Bg 

229.0 

229(20)1^ 

LAT [50] 



215.2 ±3.0 

QCD SR [51] 

Ib s / Ib 

1.236 

1.20(3)(1) 

HPQCD LAT [47] 



1.229(26) 

LAT [46] 

fo 

206.1 

204.6±5.0 

PDG [45] 

frr 

244.3 

278 ± 13 ± 10 

LAT [52] 



245(20)^2 

LAT [50] 



252.2 ±22.3 ±4 

QCD SR [53] 

fo. 

257.5 

257.5±4.6 

PDG [45] 

fo; 

272.0 

311±9 

LAT [52] 



272(16)l2 0 

LAT [50] 



305.5 ± 26.8 ± 5 

QCD SR [53] 

foj Id 

1.249 

1.258±0.038 

PDG [45] 


TABLE I: Results for the leptonic decay constants fu (in MeV). 


In the SM, the purely leptonic decays B —> £ proceed via the annihilation of the 
quark-pair into an off-shell W boson. The branching fraction for the leptonic decays is given 
by 

B(B~ -± l~v t ) = ^ fl\V ub \ 2 T B , (29) 

8vr V rnr B J 
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This work 

Data 

Ref. 

B —»• e v e 

1.16 • 10 -11 

< 9.8 • 10~ 7 

PGD [45] 



(0.88 ±0.12) • 10" 11 

UTfit [4] 



(0.85 ± 0.27) • 10" 11 

CKMhtter [54j 

B~ ->• 

0.49 • 10" 6 

< 1.0 -10” 6 

PGD [45] 



(0.38 ± 0.05) • 10” 6 

UTfit [4] 



(0.37 ±0.02) • 10” 6 

CKMhtter [54] 

B~ —>• t~v t 

1.10 • 10 -4 

(1.14 ±0.27) • 10” 4 

PGD [45] 


TABLE II: Leptonic B-decay branching fractions. 

where Gp is the Fermi coupling constant, and rri£ are the B meson and lepton masses, 
respectively, and Tp is the B meson life time. The expected branching fractions are O(10” 4 ), 
O(10” 7 ), and O(10” n ) for l — r, //, and e, respectively. The different lepton masses affect 
the values of the branching fractions through the hclicity flip factor (1 — mj/m%) 2 . 

IV. FORM FACTORS OF SEMILEPTONIC B-MESON DECAYS 

The invariant matrix element of the semileptonic decays B —> D^*'£~V£ can be written as 

M(B -A D^£~P e ) = -%V cb < D w \cO fM b\B> £0^ e , (30) 

\/2 

where the matrix elements of the semileptonic B transitions in the covariant quark 

model are defined by the diagram in Fig. [3] and are written as 
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= 7 M(1 _ 7s ) 


B(pi) 

» 



9l k + pi 


93 


fe + P2 92 


fe 


93 


0_dWo> 2 ) 


- (fc + ^13 Pi) 2 ) ^£>(*)( - ( fc + W 23P2) 2 ) 


FIG. 3: Quark model diagram for B meson semileptonic decay. 


T^ = (D(p 2 )\c0^b\B( Pl )) = 

AT ~ ~ J ^ (~{k + W 13 pi) 2 ) ( - {k + W 23 P2 ) 2 

x tr 


(2vr) 4 i " ' 

O^S 1 (k+p 1 ) 1 5 S 3 (k)'y 5 S 2 (k+p 2 ) 


= F + (q 2 ) + F_(g 2 ) q^, and 


(31) 


e\ a T^ = {D'( P2 , e2 )\cO“b\B( Pl )) = 

d 4 k 


— ^c9b9d 
x tr 


$ 


B 


{2'k)H 

0 *S 1 {k+p 1 )'f i S 3 {k) / 2 f S 2 (fc+p 2 ) 


(k + W13P1) ) $£>* - (k + W 23 P2) 


A 

“2 a 


mi + m 2 


(-g» a PqA 0 (q 2 ) + P^ 1 P a A + (q 2 ) + q^ P a A_(q 2 ) + i£ fiaPc *V( y q 2 )) . (32) 


Here, P = p 3 + p 2 , q = Pi — P 2 , and e 2 is the polarization vector of the D* meson so 
that t\ ■ P 2 = 0. The particles are on their mass-shells: p 2 = m 2 = m 2 B and p\ = m\ = 
m 2 DM . Altogether there are three flavors of quarks involved in these processes. We therefore 
introduce a notation with two subscripts Wij = m q J{m qi + rn qj ) ( i,j = 1,2,3) such that 
w ij + Wjt = 1. In our case one has q\ = b, q 2 = c, and q 3 = d. 

Our numerical results for the form factors are well represented by a double-pole 
parametrization 


Hq 2 ) 


F{ 0) _ q 2 

1 — as + bs 2 ’ m 2 


(33) 
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The double-pole approximation is quite accurate. The error relative to the exact results is 
less than 1% over the entire q 2 range. For the B transition the parameters of the 

dipole approximation are given by 



F + 

F_ 

Ao 

A + 

A_ 

V 

m 

0.78 

-0.36 

1.62 

0.67 

-0.77 

0.77 

a 

0.74 

0.76 

0.34 

0.87 

0.89 

0.90 

b 

0.038 

0.046 

-0.16 

0.057 

0.070 

0.075. 


(34) 


Since b/a is quite small for the form factors F +1 F_, A +1 A_, and V, these form factors 
show a monopole-like falloff behavior whereas A 0 has a substantial (q 2 )~ 2 contribution. In 
Fig. |4] we present our results for the semileptonic form factors within the full range of 
momentum transfer 0 < q 2 < q ^ aY , where q = (mg — m D (*)) 2 . The results of the exact 
calculations are shown by solid lines whereas the results obtained in the heavy quark limit 
are shown by dashed lines. We will discuss the heavy quark limit in the next section. It 
is interesting to note that the QCD counting rules prescribe a ( q 2 )~ l and a ( q 2 )~ 2 falloff 
behavior for the form factors F + , F_, A 0 and A + , A_, V, respectively. 

As recently noticed in 55], the ratio Fo(q 2 )/F + (q 2 ) exhibits a linear q 2 behavior 

c / 2 n , Q 2 v /„ 2 \ MQ 2 ^ 


Fo(q 2 ) = F + (q 2 ) + f-F4q 2 ), 


= 1 — aq , 


(35) 


f+(9 2 ) 

where the slope a = 0.020(1) GeV -2 was determined precisely based on lattice values of the 
two form factors. We also plot the q 2 dependence of the ratio F 0 (q 2 )/F + (q 2 ) in Fig. 0 which 
shows a linear behavior as mentioned. Our value for the slope is a = 0.019 GeV -2 which 
very well agrees with the lattice result. 
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FIG. 4: Form factors of the decays B —>• D^iv. The solid lines are the results of exact calculations 
in our approach, the dashed lines are the form factors obtained in the heavy quark limit. 
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FIG. 5: Ratio F 0 (q 2 )/F + (q 2 ). 


V. HEAVY QUARK LIMIT 


It is instructive to explore the heavy quark limit (HQL) in the heavy-to-heavy transition 
B —y D{D*). In the HQL one takes the limit m b = rn b + E, mj —> oo and rrifj = mp* = 
m c + E, m c -A oo in the expressions for the coupling constants and form factors. In this 
limit the heavy quark propagators are reduced to the static form 


S b (k+p i) = 
S c (k + p 2 ) = 


m b - ft- fa 
1 

m r 


-A 


-A 


i+ A 

+ °i 

( 1 

—2 kvi — 2 E 

\m b 

+ 

to 

+ °i 

( 1 

—2kv2 — 2 E 

\m c 


(36) 


ft~ft 2 

where pi and v t = Pi/rrii (i — 1,2) are the momenta and the four-velocities of the initial and 
final states. Moreover, we have to keep the size parameters of heavy hadrons equal to each 
other in order to provide the correct normalization of the Isgur-Wise function at zero recoil. 


By using technique developed in our previous papers, see, for instance, [56|, |57j, one can 


arrive at the following expressions for the scmileptonic heavy-to-heavy transitions defined 
by Eqs. (l3Tf and 


rpfl 

1 HQL 


t 

fc 2 v 1 HQL 


f M • z tr 1+ A)l* • 7 5 (1+ fa) = £(w) ■ « + v%), 


■ i tr c+(i+ a)j 5 - 4(1+ fa) 

• eU-g^il + w) + « + vffl 


(37) 


(38) 
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Here, w = V 1 V 2 , and the Isgur-Wise function is equal to 
J 3 (E,w) 


= 


1 00 

dr 


J 3 (E,w)= / — du$\z) [a s {z) + — a v (z) 


W 


Js(E,l )’ 

0 0 

where W — 1 + 2r(l — t)(iu — 1), z = u — 2 E^Ju/W, and 
$(z) = exp(-z/A 2 ), 





W 


(39) 


m„ 


= 


CT V (z) = 


ml + z 


ml + z ’ 

By using the definition of the form factors given by Eqs. (j3TT) and (l32i) one can easily obtain 
the expressions of the form factors in the HQL. One finds 
2 m 1 ±m 2 . . 

F±(? > = ± v^f {W> 


Ao^q 2 ) = 'JEzEt (1 + w)£(w), A + (q 2 ) =-A-(q 2 ) = V(q 2 ) = "' 1 


m 2 


ra, — m 2 . 2,ASTmJ 5(w) ' (40) 

where tc = {m\ + m| — g 2 )/(2mi7«2). We use the physical masses of the heavy hadrons 
in the numerical calculations. For the size parameter we adopt the average value A = 
(A^ + Ad + Ap*)/3 = 1.70 GeV. The parameter E characterizes the difference in mass 
between the heavy hadron and the corresponding heavy quark. We use its minimal value 
E = mr> — m c = 0.20 GeV in order to avoid the complication with confinement. 

In Fig. |4]we display the heavy-to-heavy transition form factors calculated in the HQL 
and compare them with the results of exact calculations. One can see that the two re¬ 
sults obtained with and without use of the HQL behave very similar to each other which 
demonstrates the fidelity of HQET. 

One can also consider the near zero-recoil behavior of the form factors in a similar way 
as we did in our paper on the semileptonic decay A b —* A c + tV t [40]. The standard 
parametrization of the [w — 1) expansion takes the form 


E(q\w)) = FiqU) 1 - p 2 (w - 1) + c (w 



where p 2 is called the slope parameter and c the convexity parameter. The numerical results 
are given below 



F+ 

F_ 

Ao 

^4+ 

A_ 

V 

F(qiJ 

1.12 

-0.52 

1.91 

0.99 

-1.15 

1.16 

P 2 

0.72 

0.74 

0.42 

0.93 

0.95 

0.96 

C 

0.49 

0.51 

0.28 

0.82 

0.85 

0.86 


(41) 
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which may be compared with the results obtained for the monopole form factor of a B c - 

resonance contribution: p 2 =0.71 and c = 0.51. 

It is interesting to compare the zero-recoil values of our exact form factors with the 
predictions of leading order HQET at w — 1 where £(1) = 1. One has 


E. = 


A+ = 


777-1 + 777-2 
2 y/rrii m 2 


= 1.138, 


F_ = 


7771 — 7772 

2 i y777 1 777 2 


= -0.543, 


-A. = V = + g| = 1.119, „ 0 

2^777 1 7772 777i — ?77 2 


_ 2 v /777 1 777 2 
/in — - 


= 1.993. 


(42) 


The zero-recoil values of our model form factors can be seen to be quite close to the corre¬ 
sponding HQET values except for the form factor A + where our form factor value exceeds 
the HQET result by ~ 13 %. 


VI. HELICITY AMPLITUDES AND TWO-FOLD DISTRIBUTIONS 


Let us first consider the polar angle differential decay distribution in the momentum 
transfer squared q 2 . The polar angle is defined by the angle between q = pi — p 2 and the 
three-momentum of the charged lepton k\ in the (f - ^) rest frame as shown in Fig. |HJ One 
has 


d 2 r 


|P2| V V^ni-r|2_ Gf i t ^ i2 IP 2 I 'i 7 


dq 2 d cos 6 (2tc) 3 32 




pol 


( 2*0 


IV 


cb\ 


- , 

64777? ^ 


(43) 


where |p 2 | = X 1 ^ 2 (rnl, rn 2 ,q 2 )/2mi is the momentum of the daughter meson and where we 
have introduced the velocity-type parameter v = 1 — m 2 /q 2 as well as the contraction of 
hadron and lepton tensors H IJV L nu . 


As discussed in some detail in 


40] the covariant contraction can be converted to a 


sum of bilinear products of hadronic and leptonic hclicity amplitudes using the completeness 
relation for the polarization four-vectors of the process. A synopsis of the necessary steps 
in this transformation is provided in the Appendix. 

One needs to relate the mesonic hclicity amplitudes to the invariant form factors defined 
in Eqs. (j3TT) and (J32]). To do so one requires explicit representations of the polarization 
four-vectors e^(Aiy)- They read 


= ~^(Qo,0,0, | p 2 |), e^(±) =-^=(0 ,Tl,-bO), ^(O) = -j=( 


IP 2 I ? 0, 0, qo). (44) 
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FIG. 6: Definition of angles 9, 9*, and % in the cascade decay B° —> D* + (—)• D°n + )£ v#. 

The linear relations between the two sets of form factors can then be calculated in the 
following way. 


B —D transition: 

The hclicity amplitudes are defined by H\ w = e^(Xw)T^. One obtains 


H, = 




(PgF + 


q 2 F'-) 


H± = 0, 


Pn = 


2mi |p 2 | 

—J7~ 


(45) 


Note the zero-recoil relation Hq = 0. At the other end of the spectrum at maximal recoil 
q 2 = 0 one has H t = Hq. In the Appendix we describe how to obtain the differential 
(q 2 ,cos9) distribution. One has 

dT(B De-u t ) = G 2 F \V cb \ 2 \p 2 \q 2 v 2 
dq 2 d cos 9 32(27r) 3 m^ 

x{2 sin 2 OT-Ll + 2 5e(2 cos 2 9 U L + 2 U s - 4 cos 9U SL ) }, (46) 


where we have introduced the helicity flip penalty factor 5e = m 2 /2q 2 and the hclicity 
structure functions Kl = \Hq\ 2 , dis = \H t \ 2 , and 7-Lsl = R ,e(H 0 Hj). 
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B —> D* transition: 

The helicity amplitudes are defined by H\ W \ D „ = A tJ ‘{\w)^ 2 t {^D*)T^ a - In addition to the 
iToff—shell polarization four-vectors e M (Avr) one needs the polarization four-vectors ef (A d*) of 
the D*. They read (E 2 = mi — g 0 ) 

e“(±) = -Uo , ±1, -i, 0), e“(0) = — (| P2 |, 0 , 0 , -E 2 ). (47) 

V2 m 2 

One obtains 


H t0 = e^(t)e¥(0)T^ Q = 


mi |p 2 | 


#±i±i = e^(±)el a (±)T m = 

H 00 = e^(0)eI“(0)T MQ 
1 1 


m 1 +m 2 m 2 \/q I 
1 


(P q (— A 0 + A + ) + q 2 A_) , 


m i + m 2 


P q Aq ±2m l |p 2 | V ), 


{-Pq (m\ - ml - q 2 ) A 0 + 4m] |p 2 | 2 A + ) . 


mi +m 2 2m 2 sfq* 

Note the zero-recoil relations H t0 = 0 and H±i±i = H 00 . At maximal recoil q 2 = 
dominating helicity amplitudes are H t o and H 00 with H t o = H 00 . 

The differential ( q 2 ,cos9 ) distribution finally reads (see the Appendix) 


dT{B -> D*e~v t ) G 2 F \V cb \ 2 \p 2 \q 2 v 2 


dq 2 d(cos6) 


■Br(P* Dn) 


32(2tt ) 3 m 2 

x |(1 + cos 2 6) P v + 2 sin 2 0P L — 2 cos 6 Up 
+ 2 (sin 2 6 Bp + 2 cos 2 9Pl + 2 Ps — 4 cos 6 Psl) }. 


(48) 
0 the 


(49) 


We have used the zero width approximation for the D* intermediate state which brings in 
the branching fraction Br(P* —> Dtt). The relevant bilinear combinations of the helicity 
amplitudes are defined in Table UTT1 We have dropped a factor of “3” in the definition of Ids 
and Pis compared to our paper [581] - Note that the helicity structure functions satisfy the 
zero-recoil relations 2P V — P L — P T = Pi and Pp = Pa = Ps = Psa = Pst = Ps = 0 . 
Similar relations hold for the imaginary parts. At maximal recoil one has Pp = Ps — Psl 
for the dominating helicity structure functions. 

Let us begin discussing the cos# distribution for the B -A- D*i~vn case. The distribu¬ 
tion (1461) is described by a tilted parabola whose normalized form reads 

b cos 6 + c cos 2 6 


W{6) = 


2 (a + c/3) 


(50) 
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TABLE III: Definition of helicity structure functions and their parity properties for the case 
B -> D*£~v e . 


parity-conserving (p.c.) parity-violating (p.v.) 


Hu = |ft+i+i| 2 + |ft-i-i| 2 

n L = \h 00 \ 2 

U T = Re (ft+i+itfl^) 

n IT = Im 

ft/ = 3 Re (iZ+1+itfJo + E-i-i^o) 
Un = \ Im (iL+i+iiLgo + -H’- 1 - 1 -H’o 0 ) 
ft* = l^otl 2 

ft ST = \ Re (tf+i+i/4 + 

ft/ST = I Im 

ftss = Re (ifootfJt) 

ft/SL = Im ^iLoo-Rot) 

fttot = ftt/ + fts + (Hu + fts + 3 fts) 


'Up = \H + 1 + 1 \ 2 -\H _ 1 - 1 \ 2 

Ha = ^Re [H +1+1 Hl 0 - 
ft/A = 5 Im (^+1+1^00 - H - 1-1^00) 
ftsA = £ Re (ft+i+ift^ - R-i-iRjt) 
ft/sA = £ Im (ft +1+1 ftj t - tf-1-,/4) 


The linear coefficient b/2(a + c/3) can be projected out by defining a forward-backward 


asymmetry given by 


65] 


Afb(q 2 ) — 


rfT(F) — dT( B) Jq dcos 9 dT/dcos 6 — dcos 6 1 dT/dcos 6 


dT(F) + dT(B) 


f 0 dcos 9 dr/dcos 9 + J’_ 1 dcos 9 dT/ dcos 9 
b 3 H p +AS £ Hsl 


2 (a + c/3) 


ft. 


(51) 


tot 


In the r mode there are two sources of the parity-odd forward-backward asymmetry, namely, 
a purely parity-violating source from the VA interaction leading to the ft p contribution, and 
a parity-conserving source from the VV and AA interactions leading to the %sl contribu¬ 
tion. The parity-conserving parity-odd contribution Hsl arises from the interference of the 
(0 + ; 1 _ ) and (0”; 1 + ) components of the VV and AA product of currents, respectively. In 
the case of the B —» D transition the forward-backward asymmetry arises solely from the 
(0 + ; 1 _ ) interference term of the VV product of currents. 

The coefficient c/2(a + c/3) of the quadratic contribution is obtained by taken the second 
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derivative of W(9). Accordingly we define a convexity parameter by writing 


cW) 


d 2 W(9) 

d(cos9) 2 


c , =~(i 

a + c/3 4 V 


25 P 


Ku - 2Hi 


H 


tot 


(52) 


When calculating the q 2 averages of the forward-backward asymmetry and the convexity 
parameter one has to multiply the numerator and denominator of (15T1) and (1521) by the q 2 - 
dependent piece of the phase space factor in (T46l) given by C(q 2 ) = |P 2 1 Q' 2 'C’ 2 - For example, 
the mean forward-backward asymmetry can then be calculated according to 


{•Afb) 


3 Jdg 2 C(q 2 )(H P +4 8 e H SL ) 

4 / dq 2 C(q 2 )H to t 


(53) 


Finally, integrating Eq. (146|) over cos 9 one obtains 


dT(B -> D^i~u t ) 
dq 2 


G 2 F \V cb \ 2 \p 2 \q 2 v 2 

12(2tt) 3 7tiI 


Bt(D* D7r) ■ 'Htot, 


(54) 


where "Htot = Hu + Hl + Sg (%u + Hl + 3 Hs^j ■ 

The discussion of the cos 9 distribution for the B —y case proceeds in a similar 

way except that one has to drop the contributions of the helicity structure functions Hu 
and Hp. 


VII. FOUR-FOLD ANGULAR DECAY DISTRIBUTION 


The lepton-hadron correlation function reveals even more structures when one 

uses the cascade decay B° —> D* + (—» D °to analyze the polarization of the D* 
meson. The derivation of the four-fold angular decay distribution is detailed in Appendix 
A. One has 


dT(B° -)■ D* + (D°n + )£-p e ) 
dq 2 d cos 9 d(x/2n) d cos 9* 


G 2 f lV r c &| 2 |p 2 |g 2 p 2 
(27 t) 3 12m 2 


Br(D* —>• Dtt) W(9*, 9, y), 


(55) 
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where 


W(9*,9,x) — (1 + cos 2 9) sin 2 9* Hu + ^ sin 2 9 cos 2 9* Hl — 777 cos 6 sin 2 9* Hi 

32 8 16 

9 9 

— — sin 2 9 sin 2 9* cos 2 y Ht — - sin 9 sin 29* cos x Ha 
16 8 

9 9 

H-sin 29 sin 29* cos \ Hi H— sin 9 sin 29* sin x Hu 

16 8 

9 9 

-sin 29 sin 29* sin x Hi a H-sin 2 9 sin 2 9* sin 2 % Hit 

16 16 


+ Si 


r ^ cos 2 9* Hs — ^ cos 9 cos 2 9* Hsl + ^ cos 2 9 cos 2 9* Hl 

+ — sin 2 9 sin 2 9* Hu + - sin 2 9 sin 2 9* cos 2% Ht 
16 8 

9 9 

+ - sin 9 sin 29* cos xHst — - sin 29 sin 29* cos xHi 
4 8 

9 9 

-sin 9 sin 29* sin \ His a H— sin 29 sin 29* sin x Hi a 

4 8 


9 


— - sin 2 9 sin 2 9* sin 2 % Hit 


(56) 


In our quark model all helicity amplitudes are real, which implies the vanishing of all terms 
proportional to sin y and sin2y. The angular decay distribution for the remaining terms 


agrees with the results o: 


59 


polar angle 9 used in 59421] such that 9 -A- 180° — 9. 


61] when one takes into account the different definition of the 


The four-fold distribution allows one to define a number of physical observables which can 
be measured experimentally. Integrating Eq. (156|) over cos 9* and x one recovers the two¬ 
fold ( q 2 ,cos9 ) distribution of Eq. (146]) that gives rise to the lepton side forward-backward 
asymmetry parameter Afb and the convexity parameter Cp(q 2 ). Integrating Eq. (I56p over 
cos 9 and x one obtains the hadron side cos 9* distribution described by a untilted parabola 
(without a linear term). The normalized form of the cos 9* distribution reads W(6*) = 
(a' + d cos 2 9*)/2(a! + c'/3), which can again be characterized by its convexity parameter 
given by 

d 2 W(9) d 3 Hu — 2 U L + S^Hu - 2H L - 6 H s ) 


^ d(cos9*) 2 a' + c'/3 2 H tot 

We define a normalized angular decay distribution W(9*,9,x) through 

W(9*,9, X ) 


W(9*,9, X ) = 


H 


(57) 


(58) 


tot 
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The normalized angular decay distribution W(9*,9,x) obviously integrates to 1 after 
cos 9*, cos 9, and x/27r integration. 

The remaining coefficient functions "Ht(1 — 2 Si), 1 — 2 6/), and (' Ha — 2 5{7-Lst) in 

Eq. (l56l) can be projected from the three-fold angular decay distribution Eq. (l56j) by taking 
the appropriate trigonometric moments of the normalized decay distribution W(9*,9,x)- 
The trigonometric moments are defined by 

Wi = J dcoa9dcoa9*d(x/2Tr)Mi(9*,9,x)W(9*,9,x)=< M t (9*,9,x) >, (59) 

where Mi(9*,9,x) dehnes the trigonometric moment that is being taken. One finds 

Ut 


W T (q 2 ) = < cos2x >= -2S e ) 


n 


tot 


Wi(q 2 ) = < cos 9 cos 9* cos x >= 
WA(q 2 ) = < sin 9 cos 9* cos x >= 


97t 2 (1 — 2Sg) Uj 


512 

37t 1~La — 25/Hst 

16 n 


(60) 


tot 


The coefficient functions "H^l — 28i), "Hr(l — 26e), and ( 'Ha — 28(T-Lst) can also be projected 
out by taking piecewise sums and differences of different sectors of the angular phase space 


60]. 


Finally, we consider the longitudinal and transverse polarizations of the lepton where 
we consider only the angular average of the two polarization states. For the longitudinal 
polarization one obtains 


W) = 


8{Hhf — 'Hnf _ 'Hu + Hl — $i(Hu + Hl + ^T-Ls) 


(61) 


SiHhf + Hnf T-itot 

The transverse polarization can be calculated using the representation of the polarized lepton 

tensor written down in the Appendix of 40] - One obtains 

3 TTy/di 'Hp — 2 'Hsl 


p x(q) = - 


(62) 


4^2 H tot 

For the decay B —> Dd~pg one has to drop the transverse contributions Hu and T-ip in 
Eqs (EU) and (]62]1 . It is interesting to note that for this decay there exists a very simple 
relation connecting P^.(q 2 ) and Appiq 2 ) which reads 


P*(q 2 ) = 


TTy/g 2 


2 m T 


^Fs(g 2 )- 


(63) 


The polarization of the lepton depends on the frame in which it is defined. The polariza¬ 
tion components and P^ in (1611) and (1621) are calculated in the (^~u T ) rest frame. The 


corresponding polarization components in the B rest frame have been calculated in 


62]. 
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VIII. RESULTS AND DISCUSSION 


The values of the lepton and meson masses and their lifetimes are taken from Ref. 45]. 
We also adopt the following values for the CKM matrix elements \V ub \ = 0.00413 and 
\V bc \ = 0.0411. In Fig. [Tj we represent our results for the differential branching fractions 
of the decays B —>■ D^*’iv within the full range of the momentum transfer squared. For 
comparison, we also display the form factors calculated in heavy quark limit. It is readily 
seen that both forms are very close to each other. It confirms that HQET works very well 
in the leading order for b — c transitions. In what follows we will not display the curves for 
observables obtained in the HQL. 

In Fig. |H]we represent our results for the forward-backward asymmetries of the decays 
B —> D^tv within the full range of the momentum transfer squared. The forward-backward 
asymmetry for the decay B — > Dt~D t is quite large in the lower half of the q 2 spectrum 
which can be understood from the fact that Afb = —3 5{Hsl and that 3< 5 ( g 2 ) is large in the 
threshold region. It is quite interesting that the forward-backward asymmetry for the decay 
B — > D*t~v t goes through zero at q 2 = 6.25 GeV 2 . 

The branching fractions of the decays B —>- £~u, B —> D^£~D, and B —> 7iT _ z/, as 
well as the ratios of branching fractions R(D^) are presented in Tables mm and IVl The 
branching fractions B(B —> £~v), [l = e,/i), satisfy the experimental constraints and show 
good agreement with the CKMfitter results, while the branching fraction B(B —^ r~u T ) is 
consistent with experimental data, giving more constraints on NP effects that may contribute 
to the transitions. The situation is different for the semileptonic decays. The results for 
B(B -A D^t v) are slightly larger, while the results for B(B -A D^t u t ) are slightly 
smaller in comparison with experimental data. As a result, the calculated ratios R(D ( -*' ) ) are 
slightly smaller than the SM expectation, which means they deviate from the experimental 
values even more. This may imply the appearance of NP. 

Next we define the partial hclicity rates by 


dT x G 2 f \V cb \ 2 \p 2 \q 2 v 2 


n 


dq 2 (27r) 3 I2m\ '^ X ' dq 2 4 dq 2 

where X = U, L, P,... In Figs. [9] and [H] we display the q 2 dependence of the partial differen¬ 
tial rates dTu/dq 2 , dTp/dq 2 , and the total differential rate dTu+L/dq 2 for the e mode. The 
transverse rate dominates in the low recoil region while the longitudinal rate dominates in 


dT x x dT x 

— Ol 


(64) 
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FIG. 7: Differential branching fractions of the decays B —>• D^lv. The solid lines are the results 
of exact calculations in our approach, the dashed lines are the form factors obtained in the heavy 
quark limit. 

the large recoil region. The longitudinal and thereby the total rate show a step-like behavior 
near the threshold q 2 = m 2 . Figs. flOl and fl2l show the corresponding plots for the r mode 
including the partial flip rates dTu^/dq 2 and 3dTs/dq 2 . We also show the total differential 
rate dVu + L/dq 2 + dVu + L+ 3 s/dq 2 ■ The helicity flip rates are smaller than the hclicity nonflip 
rates but contribute significantly to the total rate. 

In Figs. HU m and m we display the q 2 dependence of the convexity parameters C l F 


26 

























0.0 



FIG. 8: Forward-backward asymmetries of the decays B —>• D^i, 


v. 


and Cp for the lepton and hadron sides defined in Eqs. d52|) and (l5Tjh In the B —>■ D 
case the cos 6 distribution is described by a downward open parabola which becomes much 
flatter for the r mode. We do not plot the hadron-side convexity parameter Cp(q 2 ) for the 
B —>■ D transition since it trivially reads Cp = 3 following from the definition (J57j). For the 
B —>■ D* transition the lepton-side cos 9 distribution is again described by a downward open 
parabola which becomes almost flat for the r mode. The hadron-side cos 6* distribution is 
described by an upward open parabola which does not become flat at the zero-recoil point. 
Lepton mass effects are not very pronounced. 
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Unit 

This work 

Data 

Ref. 

B° -> D+i~P 

10" 2 

2.74 (2.65) 

2.17 ±0.12 

HFAG [9] 




2.21 ± 0.16 

BABAR [63] 

B° -> D+t~u t 

10~ 2 

0.73 (0.71) 

1.02 ±0.17 

BABAR [10] 

B° -> D*+tTv 

10~ 2 

6.64(7.21) 

5.05 ±0.12 

HFAG [9] 




5.49 ± 0.30 

BABAR [63] 

B° D*+t~v t 

10' 2 

1.57(1.70) 

1.76 ±0.18 

BABAR [10] 

B° -> ir+e-jj 

10" 4 

1.69 

1.41 ±0.09 

BABAR [7] 




1.49 ±0.08 

Belle [8] 

B° — > tt + t~u t 

10~ 4 

1.01 




TABLE IV: Semileptonic decay branching fractions of B meson. The values obtained in the HQL 
are given in brackets. The experimental errors are combined in quadrature. 



This work 

SM 

Data 

R(D) 

0.265(0.268) 

0.297 ±0.017 

0.388 ±0.047 

R(D*) 

0.237 (0.235) 

0.252 ±0.003 

0.321 ± 0.021 


TABLE V: Ratios of branching fractions R(D ) and R(D*) calculated in our model (the values ob¬ 
tained in the HQL are given in brackets) and compared with the SM expectations and experimental 
data. 


In Figs, m DU and H%1 we show plots of the q 2 dependence of the longitudinal, transverse 
and total polarization of the lepton for the B —> transition. In the case of the electron 

the curves reflect the chiral limit of a massless lepton in which the lepton is purely left- 
handed, i.e. one has Pf = —1, Pf = 0, and |P f | = 1. For l = r the transverse polarization 
is large and positive and dominates the total polarization. The transverse polarization of 


the r drops out after the appropiate azimuthal averaging, as has been done in 


13]. Note 


that the transverse polarization in the r mode results solely from the scalar-longitudinal 
interference contribution BisL- The longitudinal polarization has switched its sign relative 
to the rri£ = 0 case. 














The corresponding curves for the B —y D*£ ry transition are shown in Figs. [19], [20] and 
EU The longitudinal and transverse polarization components are distinctly different from 
their rri£ = 0 values Pf = —1 and Pj = 0. The longitudinal component becomes larger 
in magnitude when q 2 increases while the transverse polarization becomes smaller as q 2 
increases. At zero recoil the transverse polarization of the charged lepton PJ tends to zero 
in agreement with the vanishing of Bp and Bsl at zero recoil. The total polarization of the 
r shown in Fig. [21] has an almost flat behavior with \P e \ ~ 0.7. The overall picture is that 
the polarization is mostly transverse at threshold and turns to longitudinal as q 2 reaches the 
zero-recoil point. 

In Figs. [221 ESI and EH we display the q 2 dependence of the three trigonometric moments 
W t (i = T, /, A) of the normalized three-fold angular function W ( 9 *, 9, y) defined in Eq. (16011 . 
Lepton mass effects can be seen to be quite large for all three moments. 

Finally, in Figs. 1251 and 1261 we present the q 2 dependence of the rate ratios (£ = e,/i) 

dF(B ->■ D^t~u t ) / dT(B D^iTVt) 

dj 2 / d,q 2 ‘ ( ’ 

Hopefully there will be enough data in the future to explore the apparent flavor violation in 
the tauonic semileptonic B —>■ D^*> transitions in more detail by measuring the rate ratios 
in different q 2 bins. 

Next we present our model results for the average values of the polarization observables: 
the forward-backward asymmetry < Aps >, the convexity parameter < Cp >, the leptonic 
< z > polarization components, and the three trigonometric moments < W t (i — T, /, A). 
Lepton mass effects can be seen to be quite large for the average values of the polarization 
observables. 


-^£> 0 ) (Q 2 ) 
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FIG. 9: B —>■ D transition: the g 2 dependence of the partial rate dr^/dg 2 for the e mode (in 
units of 1(T 15 GeV -1 ). 




FIG. 10: B —>■ D transition: the g 2 dependence of the partial nonflip rates dr^/dg 2 , and the 
flip rates dTu^/dq 2 and 3dTs/dq 2 for the r“ mode (in units of 10 -15 GeV -1 ). Also shown is 
the total rate dT^/dg 2 + dTi/dq 2 + 3 dTs /dg 2 . 
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FIG. 11: B —> D* transition: the q 2 dependence of the partial rates dTu/dq 2 (dashed), dr L/dq 2 
(dot-dashed) and their sum dTu + z J /dq 2 (solid) for the e~ mode (in units of 10~ 15 GeV -1 ). 




FIG. 12: B —> D* transition: the q 2 dependence of the partial nonflip rates dTjj^/dq 2 , and 
the flip rates dTu^/dq 2 and 3dTs/dq 2 for the t~ mode (in units of 10 -15 GeV^ 1 ). Also shown 
is the total rate dVu+L/dq 2 + dVjj+L/dq 2 + 3 dTs/dq 2 . 
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FIG. 13: B —> D transition: the q 2 dependence of the lepton convexity parameter Cp(q 2 ) for 
the e~- (solid) and r _ -mode (dashed). 



FIG. 14: B —> D* transition: the q 2 dependence of the lepton convexity parameter Cp(q 2 ) for 
the e~- (solid) and T - -mode (dashed). 
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FIG. 15: B —> D* transition: the q 2 dependence of the hadron convexity parameter Cp(q 2 ) for 
the e~- (solid) and T~-mode (dashed). 



FIG. 16: B —>■ D transition: The q 2 dependence of the longitudinal polarization component 
P~(q 2 ) for the charged leptons e~- (solid) and r _ -mode (dashed). 
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FIG. 17: B —> D transition: the q 2 dependence of the transverse polarization component P^.(q 2 ) 
for the charged leptons e~- (solid) and r~-mode (dashed). 



FIG. 18: B —>■ D transition: the q 2 dependence of the total lepton polarization \P^\{q 2 ) = 
(P%) 2 + (Pf ) 2 for the e~- (solid) and r“-mode (dashed). 
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FIG. 19: B —> D* transition: the q 2 dependence of the longitudinal polarization component 
P~{q 2 ) for the charged leptons e~- (solid) and r _ -mode (dashed). 



FIG. 20: B —> D* transition: the q 2 dependence of the transverse polarization component 
Px{q 2 ) for the charged leptons e~- (solid) and r _ -mode (dashed). 
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FIG. 21: B —> D* transition: the q 2 dependence of the total lepton polarization \P e \(q 2 ) = 
\/ (P£) 2 + (Pf) 2 for the e~- (solid) and r _ -mode (dashed). 



FIG. 22: B —>■ D* transition: the q 2 dependence of the trigonometric moment Wt dehned in 
Eq. (l60l) for the e~- (solid) and T _ -mode (dashed). 
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FIG. 23: B —> D* transition: the q 2 dependence of the trigonometric moment Wj defined in 
Eq. (1601) for the e~- (solid) and r _ -mode (dashed). 



FIG. 24: B —> D* transition: the q 2 dependence of the trigonometric moment Wa defined in 
Eq. (1601) for the e~- (solid) and T~-mode (dashed). 
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FIG. 25: The q 2 dependence of the ratio R{D) 
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FIG. 26: The q 2 dependence of the ratio R(D*) 









B - 

> D 





r L 

e 

11.9 







T 

1.05 








fV 

rv 

r sl 





T 

0.25 

0.62 

0.38 








B -> 

■ D* 





rv 

rv 

T t 

I V 

Tp 

Ta 


e 

13.2 

15.6 

5.35 

8.94 

-7.42 

-3.01 


T 

3.02 

2.08 

1.32 

1.70 

-1.42 

-0.44 



IV 

?V 

rv 

f/ 

fs 

f SL 

IVt 

T 

0.64 

0.46 

0.27 

0.37 

0.20 

0.29 

0.22 


TABLE VI: q 2 averages of the rate functions in units of 10 15 GeV. We do not display the helicity 
flip results for the e mode because they are of the order of 10 -6 — 10“' in the above units. 
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B -> D 



< Ap B > 

<c e F > 

<c h F > 

e 

-1.17 (-1.16) • 10~ 6 

-1.5 (-1.5) 

3(3) 

r 

-0.36 (-0.36) 

-0.26 (-0.26) 

3(3) 


< P z > 

< P x> 

< \ pe \ > 

e 

-l(-l) 

0(0) 

1(1) 

T 

0.33(0.33) 

0.84 (0.84) 

0.91(0.91) 


B - 

> D* 



< A ( FB > 

<c e F > 

A 

V 

e 

0.19(0.18) 

-0.47 (-0.44) 

0.93 (0.88) 

T 

0.027(0.021) 

-0.062 (-0.057) 

0.58 (0.52) 


<Pf > 

< P x > 

< \P l \ > 

e 

-l(-l) 

0(0) 

1(1) 

T 

-0.50 (-0.51) 

0.46 (0.43) 

0.71 (0.71) 


< Wt > 

< W f > 

< W A > 

e 

-0.093 (-0.098) 

0.054 (0.055) 

0.062 (0.059) 

T 

-0.057 (-0.059) 

0.025 (0.025) 

0.077 (0.074) 


TABLE VII: q 2 averages of polarization observables. For comparison with results from the HQL, 
we add in brackets the corresponding HQL values. 
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IX. SUMMARY AND CONCLUSIONS 


We have provided a detailed analysis of the pure leptonic and semileptonic decays B —> 
and B —* D^*H~ U£ (i = e, fi, t) within the SM in the framework of our covariant quark 
model with built-in quark confinement. We have described in some detail how to compute 
the one-loop quark contributions needed for the calculation of the transition form factors 
including a discussion of how the confinement of the constituent quarks is achieved in the 
covariant quark model. In the light of the recent experimental indications for a possible 
breaking of lepton universality in the r sector we have put particular emphasis on how to 
isolate heavy lepton mass effects in the semileptonic decays. 

We have described how to obtain the full angular decay distributions for B —> Dl~ U( and 
the cascade decay process B —> D*(—> Dtt)£~vh as well as the corresponding angular decay 
distributions for their charge-conjugate processes. The coefficients multiplying the angular 
factors in the angular decay distributions have been given in terms of hclicity structure 
functions for which we have provided simple expressions for the maximal recoil q 2 = 0 and 
the minimal (zero) recoil q 2 = (m i - m 2 ) 2 . Starting from the angular decay distributions we 
have defined a multitude of polarization observables for which we have provided numerical 
results on their q 2 spectra and their q 2 averages for zero and nonzero lepton masses. The 
polarization observables include the transverse and longitudinal polarizations of the charged 
t~ which considerably deviate from their simple rri( = 0 left-chiral structure. 

We are looking forward to a wealth of data on these decays expected in the near future 
which will allow one to deeply probe into their decay structure, in particular for the tauonic 
mode. Such an analysis will reveal possible deviations from the SM predictions not only in 
the branching fractions of the processes but also in the multitude of polarization observables 
and their q 2 spectra. 
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Appendix A: Spin kinematics 


In this Appendix we provide a synopsis of how to obtain the angular decay distributions 
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for the decays B —» D^*>i zy following the description in 
tion of the angular decay distribution is given by 

W'(0‘,6, X ) = hep Pf‘Pf 0 PZlPm L„ v , 


64,]. The covariant representa- 


(Al) 


where h a >p> is the hadronic decay tensor for the decay D* —* D + 7 r (with h a P = P 3 P 3 /1P 3 1 2 ), 
H afl H\ v is the tensor describing the decay B —» D* + W~^_ shell , and is the lepton tensor 
describing the decay W' 0 ^_ shell —> £~ + &£. The tensors are connected by the appropriate spin 
1 and spin ( 0 © 1 ) pro pag ator projectors P^" (q) = —g^ + q^q^/q 2 and P^ 1 (g) which, in the 
unitary gauge, reads [ 66 ] 


M)©! 


( 9 ) = - g 1,0 + 


■jV = f _ < ei + ‘Pi 0 


m 


w 


q‘ 


_ q 2 
q 2 


m 


w 


—^— 

spin 1 


spin 0 


= Pf(q) - (1 - \)Pf(q), 
m w 


(A 2 ) 


(A3) 


where PQ U (q) = q^q v /q 2 is the spin 0 propagator. The factor (1 — q 2 /rriyy) multiplying the 
spin 0 propagator in Eq. 0A2j) is usually set to 1 in low energy applications as in the decay 
B —y D*£~U£. For example, at the highest q 2 value at q 2 = (m b — mp *) 2 the correction 
amounts to a mere 0.17% and will therefore be dropped in the following. 

In order to convert the covariant representation of the angular decay distribution Eq. (ED 
to the hclicity representation one makes use of the completeness relations 


PZ 1 ( 9 ) = -sT + 


qVqP 


m 


— ^ £ tJ ‘(rn)£^ u (m')g ri 

W m,m'=t,±, 0 

Pi V (q) = —g^ u + q^q u /q 2 = ^ £^{m)£ fv {m'), 

m,m'= d=,0 


(A4) 

(AS) 


where the tensor g mm > = diag(+, —, —, —) is the spherical representation of the metric tensor 
whose components are ordered in the sequence m, m! = t, +, 0, —. With the help of the com¬ 
pleteness relation (lASjl one can convert the covariant form of the angular decay distribution 
Eq. (lAip into the hclicity form 

(-iy +J 'Sx 2 , w s>,x w KM{e‘)H X 2 Kw (j)Hi A , f,e,x). 


w\e\e, x )= Y 


(A 6 ) 
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In (1A6D we have chosen a representation of the hclicity amplitudes which is particularly well 
suited for computer processing. Compared to the helicity amplitudes introduced in the main 
text we have used H 0 Xw =o (J = 0) = H ot and H 0i±1 x w =o,±( J = 1) = #o,±i 0 ,±i. 

The hclicity representation of the hadronic decay tensor h a p(9*) describing the decay 
D* —> Dtt is given by 


/ 


hx 2 x'(0*)=d 1 OX2 (9*)d 1 ox ,(e*) = 


\ sin 2 9* 


^ sin 20* 


A sin 2 9* 


\ 


+ 


2V2 


sin 29* 


cos 2 9* — 273 sin 29* 


V 


- \ sin 2 9* sin 29* \ sin 2 9* 




( 


(1 =F cos 9f 


Sill 


2 9e 2ix ^ 


T 


2 sin 2 9 


V 


^=(1 =F cos 9) sin0e * x 

sin 2 9e~ 2lx =f^=( 1 ± cos#) sin0e 

( o„;„2 a 2 o a„iy o \ 


3—®x 


^(1 ± cos0) sin 9e %x 
(1 ± cos6*) 2 


2 sin 2 9 
— sin 2 9e~ lx 


- sin 2 9e lx 
4 cos 2 9 


V 


-2 sin 2 9e~ 2ix sin 29e~ ix 


-2 sin 2 9e 2ix 
sin 2 9e lx 
2 sin 2 9 


(AT) 


For the hclicity representation of the lepton tensor one obtains (v — l — m 2 /q 2 ) 

(2q 2 v)- 1 L Xw y w (l,l,9,x) = 

=F^(1 =F cos 9) sin 9e ix 


40] 


(AS) 


The upper/lower signs in the nonflip part of (1A8|) stand for the (£ vg) case relevant for the 
decays B° ->■ D^+£-pg and B~ D^°t ug, and the (£ + ui) case relevant for the decays 
B+ D^ 0 £+u e and B° —> £ + ug . The spin 0/spin 1 interference contribution is given 

by 


(2g 2 u) 1 L OXw {0,1 , 9, x) = 


( 2 <? v) L* Xw 0 (l, 0, 9, x) = 


St. ( — ^sin#e lx 4 cos 9 ^ sin 9e lx j , (A9) 


(A w = 1, 0, — 1) and 


(2q 2 v)~ 1 L oo (0, 0, 9, x) — 4 8 e . (A10) 

For the cos 9 distribution of the decay B —> D*i~ug written down in Eq. (I40H one needs the 
integrated form of Eq. (1A6D . One obtains 


W\9) = j dcos 9*d X /2nW'(9*, 9, X ) = 

2 V (-1 ) J+r Sx 7 x w H X7Xw (J)H- Mw (J')L XwXw (J,J\9), (All) 

J, J',Xw^2 


43 












where 


L\ W \ W (J, J\ 0) = I dx/2irL Xw \ w (J, J', 9 , y). (A12) 

The integration (IA12D is easily done. The result is given by Eqs. (1A81A9D where all terms 
proportional to e ±lx ,e~ ±2lx have been dropped. The cos 9 distribution for B —> D£~D^ 
written down in Eq. (|46li is obtained from (1A1 1 j) by omitting S\ 2 \ w and dropping the label A 2 . 
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